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Abstract 



Q^ , We consider the wireless two-way relay channel, in which two-way data transfer takes place between 

the end nodes with the help of a relay. For the Denoise-And-Forward (DNF) protocol, it was shown by 
Koike-Akino et. al. that adaptively changing the network coding map used at the relay greatly reduces the 
impact of Multiple Access interference at the relay. The harmful effect of the deep channel fade conditions 
can be effectively mitigated by proper choice of these network coding maps at the relay. Alternatively, 
^^ ■ in this paper we propose a Distributed Space Time Coding (DSTC) scheme, which effectively removes 

^ ' most of the deep fade channel conditions at the transmitting nodes itself without any CSIT and without 

CN ' any need to adaptively change the network coding map used at the relay. It is shown that the deep 

^f) , fades occur when the channel fade coefficient vector falls in a finite number of vector subspaces of 

CO . C^, which are referred to as the singular fade subspaces. DSTC design criterion referred to as the 

CS| ' singularity minimization criterion under which the number of such vector subspaces are minimized is 

obtained. Also, a criterion to maximize the coding gain of the DSTC is obtained. Explicit low decoding 

. . , complexity DSTC designs which satisfy the singularity minimization criterion and maximize the coding 

rS 

^ . gain for QAM and PSK signal sets are provided. Simulation results show that at high Signal to Noise 

C^ ■ 

Ratio, the DSTC scheme provides large gains when compared to the conventional Exclusive OR network 

code and performs slightly better than the adaptive network coding scheme proposed by Koike-Akino 

et. al. 

I. Background and Preliminaries 

A. Background 

We consider the two-way wireless relaying scenario shown in Fig{T] Two-way data transfer takes place 
between the nodes A and B with the help of the relay R. It is assumed that all the three nodes operate in 
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half-duplex mode, i.e., they cannot transmit and receive simultaneously in the same frequency band. The 
idea of physical layer network coding for the two way relay channel was first introduced in |[T], where the 
multiple access interference occumng at the relay was exploited so that the communication between the 
end nodes can be done using a two phase protocol. A protocol called Denoise-And-Forward (DNF) was 
proposed in ||2l, which consists of the following two phases: the multiple access (MA) phase (Fig. 1 1(a) i. 



during which A and B simultaneously transmit to R and the broadcast (BC) phase (Fig. |l(b)[ ) during 
which R transmits to A and B. Network coding map, which is also referred to as the denoising map, is 
chosen at R in such a way that A (B) can decode the messages of B (A), given that A (B) knows its own 
messages. During the MA phase, the transmissions from the end nodes were allowed to interfere at R, 
but the harmful effect of this interference was mitigated by a proper choice of the network coding map 
used at R. Information theoretic studies for the physical layer network coding scenario were reported in 
im, H. A differential modulation scheme with analog network coding for bi-directional relaying was 
proposed in O. The design principles governing the choice of modulation schemes to be used at the 
nodes for uncoded transmission were studied in |!6l. An extension for the case when the nodes use 
convolutional codes was done in [7|. A multi-level coding scheme for the two-way relaying scenario was 
proposed in ||8]. Power allocation strategies and lattice based coding schemes for bi-directional relaying 
were proposed in ||9]- 

It was observed in ||6l that the network coding map used at the relay needs to be changed adaptively 
according to the channel fade coefficients, in order to minimize the impact of the Multiple Access Inter- 
ference (MAI). A computer search algorithm called the Closest-Neighbour Clustering (CNC) algorithm 
was proposed in @ to obtain the adaptive network coding maps resulting in the best distance profile 
at R. An adaptive network coding scheme for MIMO two-way relaying based on the CNC algorithm 
was proposed in |[TOl . An alternative procedure to obtain the adaptive network coding maps, based on 
the removal of deep channel fade conditions using Latin Squares was proposed in iTTl. A quantization 
of the set of all possible channel realizations based on the network code used was obtained analytically 
in |[T2l . An extension of the adaptive network coding scheme for MIMO two-way relaying using Latin 
Rectangles was made in |[T3l . 

As an alternative to the adaptive network coding schemes in ||6l and ifTTI - |[T2l . in this paper, we 
propose a Distributed Space Time Coding (DSTC) scheme, which mitigates the effect of MAI to the 
fullest extent possible at the transmitting nodes itself without any CSIT For the proposed DSTC scheme 
the network coding map used at R need not be changed adaptively according to channel conditions which 
reduces the complexity at R to a great extent and also eliminates the need for overhead bits from R to 

February 26, 2013 DRAFT 



A and B to indicate the choice of the network coding map. 

A distributed space time coding scheme for a wireless two-way relay network with multiple relay 
nodes was proposed in |[T4l . in which the DSTC was constructed at the relay nodes. In the proposed 
scheme, the DSTC is constructed at the end nodes A and B. 

B. Signal Model 

Throughout, a quasi-static fading scenario is assumed with the Channel State Information (CSI) 
available only at the receivers. Let Iia and Kb denote the fade coefficients associated with A-R and 
B-R links and /i'^ and h'^ denote the fade coefficients associated with R-A and R-B links. All the fading 
coefficients are assumed to follow Rician distribution. 

Let S denote the unit energy M = 2^ point constellation used at the end nodes. Let /x : Fg — >• 5 
denote the mapping from bits to complex symbols used at A and B. 

1) Denoise-And-Forward (DNF) protocol: In the sequel, we briefly describe the adaptive network 
coding schemes based on the DNF protocol proposed in ||6], lOTI - |[T2l . Throughout the paper, by DNF 
protocol, we refer to the schemes proposed in |0 and lITTI - [|T2l . 

In the DNF protocol, transmission occurs in two phases: Multiple Access (MA) phase during which 
A and B simultaneously transmit to R and Broadcast (BC) phase during which R transmits to A and B. 

MA Phase: Let xa = fJ-{sA), xb = IJ'{sb) G S denote the complex symbols transmitted by A and B 
respectively, where sa,sb G Fg. The received signal at R is given by, 

yn = hAXA + fiBXB + ZR. 

The additive noise zr is assumed to be CM{0,a^), where CJ\f{0,a'^) denotes the circularly symmetric 
complex Gaussian random variable with mean zero and variance a^. 

BC Phase: Let {xa, xb) G 5^ denote the Maximum Likelihood (ML) estimate of {xa, xb) at R based 
on the received complex number y^. Depending on the value of Ha and Kb, R chooses a many-to-one 
map Jvi^A'^B . ^2 _i. ^/^ where S' is the signal set (of size between M and A'P) used by R during the 
BC phase. 

In order to ensure that A (B) is able to decode B's (As) message, the map _M^-*'^b should satisfy the 
exclusive law |[6||, i.e., 

M^^'^''{xa,xb) ^ M''^'^''{x'^,Xb), foTXAJ^x'^, \/xB^S, 

M'^^'''''{xa,xb) ^ M''^''^''{xa,x'q), ioiXB^x'j^, V x^ G 5. 
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The CNC algorithm proposed in [[6] obtains the map Jvl^^^hs ^vhich results in the best distance profile 
during the MA phase at R. The CNC algorithm is run for all possible channel realizations and a partition 
of the set of all channel realizations is obtained depending on the chosen network coding map. For a 
given channel realization, the choice of the network coding map is indicated to A and B using overhead 
bits. During the BC phase R transmits xr = M^^'^'^{xa,xb) G S' . The received signals at A and B 
during the BC phase are respectively given by, 

VA = h'j^XR + ZA-, 
VB = h'BXR + ZB, 

where za and zb are independent and CM{0,a^). Since the map Jli^'--^'^^ satisfies the exclusive law and 
A (B) knows its own message xa (xb), it can decode xb (xa) by decoding xr. 

The CNC algorithm optimizes the entire distance profile instead of maximizing only the minimum 
distance. In some cases, this results in the use of signal sets with a larger cardinality during the BC 
phase. To solve this problem, an algorithm called the Nearest Neighbour Clustering (NNC) algorithm 
was proposed in ||6l which maximizes the minimum distance alone, instead of optimizing the entire 
distance profile. 

The choice of the network coding map obtained depends only on the ratio ^ and not the individual 
values of Ha and Hb I©. In ifTTI . the values of ^ for which deep channel conditions occur were 
identified and network coding maps which remove the harmful effect of these deep channel conditions 
were obtained by the completion of partially filled Latin Squares. 

2) The Proposed DSTC Scheme: For the proposed DSTC scheme, transmission occurs in four phases: 
Two MA phases during which A and B simultaneously transmit to R followed by two BC phases during 
which R transmits to A and B. Two independent complex symbols each from A to B and B to A get 
exchanged at the end of the four phases and hence the information rate in bits per channel use for the 
proposed scheme is same as that of the DNF protocol. 

MA Phases: Let xa^ = /u(s^J, x^^ = ^(^^2) ^ "^ denote two independent complex symbols A wants 
to communicate to B. Similarly, B wants to communicate two independent complex symbols xb^ = 
fi{sBi),XB2 = /^(•Sfia) G 5 to A. During the i^^ MA phase i G {1,2}, A transmits /^(x^^jX^J G C, a 
function of xa^ and x^^i ^^^ similarly B transmits /]j (x^^ , x^^ ) G C, a function of xb^ and x^^. The 
received signal at R during the two MA phases can be written as, 
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YR = [2/R1 y/jj = [hA hs] 



fhixB,,XB2) fU^B,,XB2) 



+ 



[ ZRi 



ZR2 



where i/r. denotes the received signal at R during the i*^ MA phase, zr^ and zr^ are independent and 
C7V(0, (T^). Let xa = [xAja^^a] ^^i^ ^B = [^^SiX^J. The matrix, 



C(xa,xb) 



(1) 



represents a DSTC codeword matrix. Note that in the DSTC codeword matrix, x^^ and x^^ can occur 
only in the first row and, xri and XB2 can occur only in the second row. In this way the DSTC differs 
from space time codes for the conventional 2x1 multiple antenna system with two collocated antennas 
at the transmitter in which the complex symbols can occupy any entry in the codeword matrix. 

For a complex number x, let x^ and x^ denote the real and imaginary parts of x. 

Definition 1: A DSTC is said to be linear if the entries of the first row of the codeword matrices are 
complex linear combinations of x^ , x^ , x^ , x\ and the entries of the second row are complex linear 
combinations of x^ ,x|j ,x|j ,x|j . Any codeword matrix C(xa,xb) of a linear DSTC can be written 
as. 



C(xa,xb) = Y. ^lA. + ^\A. + Wg^^g, + ^"b/b,- (2) 

i=l,2 

The matrices W^. , W^.,Wg. and Wg. are referred to as the weight matrices of the DSTC. Note 
that the entries of the second (first) row are zeros in the matrices W^. and W^ (Wg and Wg ). 

Definition 2: A linear DSTC is said to be over the signal set S if the entries of the first (second) row 
of the codeword matrices are complex linear combinations of xa-, and XA2 (xbi and XB2), where xa^, 
XA2) xbi and XB2 belong to the signal set S. 

xaMa 
For a linear DSTC over S, codeword matrix C(xa,xb) is of the form C(xa,xb) = 

[xbMb 

where Ma and Mb are 2 x 2 complex matrices referred to as the generator matrices at node A and B 

respectively. Throughout the paper, we consider only linear DSTCs over a signal set S. 

BC Phases: Let {sai,sa2,sbi,sb2) denote the maximum likelihood estimate of {sai,sa2,sbt_,sb2) at 

R. The relay R transmits xr^ = ^{sai © sbJ and xr^ = fJ-{sA2 © SB2) during the first and second BC 

phases respectively, where © denotes the bit-wise XOR operation. The received signals at the end nodes 

during the two BC phases are given by, yA, = h'j^XR. + za, and y^, = h'gXR^ + 2;^, , where i G {1, 2}. 
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Since A (B) knows its own messages and the XOR map satisfies the exclusive law, A (B) can decode 
SB, (saJ i G {1,2}, by decoding xr^. 

Note that for the proposed DSTC scheme the signal set used during the BC phase is of the minimum 
cardinality 2^ (the cardinality of the signal set should be at least 2'*' to convey A information bits). In 
contrast, for the scheme proposed in ||6], depending on channel conditions unconventional signal sets 
with cardinality greater than the minimum cardinality are required. Minimum cardinality signal set is 
used during the BC phase and throughout the paper the focus is on optimizing the performance during 
the MA phase. 

Some of the advantages of the proposed DSTC scheme over the schemes proposed in ||6], lITTI - |[T2l 
are summarized below: 

• Unlike the schemes proposed in ||6], lITTI - |[T2l . for the proposed DSTC scheme, the network coding 
map used at R need not be changed adaptively according to channel conditions. Any network coding 
map satisfying the exclusive law will give the same performance and for simplicity, the conventional 
bit-wise Exclusive OR (XOR) map itself can be used. This reduces the complexity at R to a great 
extent and also eliminates the need for overhead bits from R to A and B to indicate the choice of 
the network coding map. 

• For the scheme proposed in ||6], for certain channel conditions the adaptive network coding map 
necessitates the use of unconventional signal sets with cardinality greater than the minimum cardi- 
nality required during the BC phase, which results in a degradation in performance. For the proposed 
scheme, the relay always uses a conventional signal set with minimum cardinality. 

• The adaptive network coding maps were obtained in |[6l, by exhaustive computer search. For 
the proposed scheme no such computer search is required, since the same network code is used 
irrespective of channel conditions. 

The contributions and organization of the paper are as follows: 

• For a classical rit x n,. MIMO system with collocated antennas, deep channel fade conditions occur 
when the channel fade coefficient vector belongs to a finite number of vector subspaces of C"* 
referred to as the singular fade subspaces. The way in which transmit diversity schemes (space 
time codes) remove the harmful effect of these singular fade subspaces is discussed. The connection 
between the dimension of these singular fade subspaces and the transmit diversity order is explained 
(Section II). 

• The MAC phase of the DNF protocol for the two-way relaying scenario can be viewed as a virtual 
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2x1 MISO system. The singular fade subspaces for the classical 2x1 MISO system, are singular 
fade subspaces for the two-way relaying scenario as well. The connection between dimension of these 
singular fade subspaces and the diversity order for the adaptive network coding schemes proposed 
in la and O- El is discussed (Section III A). 

• The singular fade subspaces for the proposed DSTC scheme are identified. The goal of minimizing 
the number of singular fade subspaces results in a new design criterion referred as the singularity 
minimization criterion for DSTCs. It is shown that for a properly chosen DSTC, most of the vector 
subspaces which were singular fade subspaces for the DNF protocol, are no longer singular fade 
subspaces for the DSTC scheme. Also, a criterion to maximize the coding gain of the proposed 
DSTC scheme is obtained (Section III B). 

• It is shown that for DSTCs which are over S, where 5 is a square QAM or 2'^-PSK signal set, 
the coding gain is maximized when the generator matrices Ma and Mb at nodes A and B are 
unitary matrices. Explicit construction of DSTCs over QAM and PSK signal sets which satisfy the 
singularity minimization criterion and maximize the coding gain are provided. It is shown that for 
all DSTCs over S with unitary generator matrices Ma and Mb, the ML decoding complexity at 
R is 0{M^) for any arbitrary signal set and is 0{M'^) for square QAM signal sets. Note that the 
brute force ML decoding complexity is 0{M^) (Section IV). 

• Simulation results presented in Section V show that at high SNR, the DSTC scheme provides large 
gains when compared to the conventional XOR network code based on the DNF protocol and 
performs slightly better than the adaptive network coding scheme proposed in ||6l. 

Notations: The complex number \/^ is denoted by j. The set of integers, Gaussian integers, rational, 
real and complex numbers are respectively denoted as Z, Z[j],Q, R and C All the vector spaces and 
vector subspaces considered in this paper are over the complex field C, unless explicitly mentioned 
otherwise. Throughout, vectors are denoted by bold lower case letters and matrices are denoted by bold 
capital letters. Let CA/'(0, a^In) denote the circularly symmetric complex Gaussian random vector with 
zero mean and covariance matrix cr^In, where In denotes the n x n identity matrix. Let (ci, C2, . . . cl) 
denote the vector subspace over C spanned by the complex vectors ci, C2, . . . cl- For a matrix A, A^ 
and A^ denotes its transpose and conjugate transpose respectively. For a vector subspace F of a vector 
space, V-^ denotes the vector subspace {x : x^ v = 0, Vv G V} and dim(V^) denotes the dimension of V. 
The all zero vector of length n is denoted by On- For a square matrix A, let rank(A) denote its rank and 
let det(A) denote its determinant. For a complex number x, x^ and x^ denote the real and imaginary 
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parts of X, X* denotes its conjugate and |x| denotes its absolute value. For a vector v, || v || denotes its 
Euclidean norm. For a matrix A, Row{A) and Col{A) respectively denote the row space and column 
space of A. K{X) denotes the expectation of X. 

II. The Notion of Singular Fade Subspaces for the Collocated MIMO system 

In this section, to explain the notion of singular fade subspaces, we digress from the two-way relaying 
scenario and focus on the classical MIMO system with collocated antennas. Consider the classical MIMO 
system with nt transmit antennas at the transmitter Tx and n^ receive antennas at the receiver Rx, with 
H being the n^ x nt complex fade coefficient matrix. The entries of the matrix H are assumed to be 
i.i.d. and Rician distributed. 

A. Singular Fade Subspaces for the Collocated MIMO system with Spatial Multiplexing 

Consider the spatial multiplexing of independent complex symbols at Tx, i.e., the received complex 
vector at Rx is given by y = Hx + z, where x is the transmitted message vector of length nt whose 
components independently take values from the signal set S and z is CAA(0, cr^lnj. 

Let cSr^(H) C C"-- denote the effective signal set at Rx, i.e., 5fe(H) = {Hx : x G S"^*}. Let 
A5 denote the difference constellation of the signal set S, i.e., A5 = {s — s' : s,s' G S]. The 
distances between two points in the effective constellation Srx{^) are of the form || HAx ||, where 
Ax/On^,AxG A5"*. 

Definition 3: For an nt x n^ MIMO system, the channel fade coefficient matrix H is said to be a 
deep fade matrix if the minimum distance of the effective constellation Srx (H) is zero. The row space 
of a deep fade matrix is said to be a deep fade space. 

Let hk, ^ <k < nr, denote the k^^ row of H. Since || HAx p= X]fe=i l^k^xp, for the minimum 
distance of the effective constellation Sjix(H.) to be zero, all the vectors hj, 1 < fc < n^, should fall in 
a vector subspace of the form (Ax)-*- for some Ax G A5"*. In other words, for || HAx || to be zero, 
the row space of H should be a subspace of the vector subspace of C"* of the form (Ax)-*- for some 
Ax G A5"*. The vector subspaces of the form (Ax)^ are referred to as the singular fade subspaces. 
Formally, a singular fade subspace can be defined as follows: 

Definition 4: A vector subspace V of C"' is said to be a singular fade subspace if all the vector 

subspaces of V are deep fade spaces. 





Note that 



is always a singular fade subspace referred to as the trivial singular fade subspace. 



February 26, 2013 DRAFT 



Example 1: Consider the 2x1 MISO system with spatial multiplexing with 4-PSK signal set S = 
{±1, ±j}. The difference constellation of 4 PSK signal set has 9 points AS = {0, ±2, ±2j, ±1 it j}. For 
this case, the set of fourteen singular fade subspaces, which are of the form (Ax)-*-, where Ax G A5^ 
are given by, 
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(3) 

The fade coefficient matrix (which is a row vector for this example) is a deep fade matrix (vector) if the 

row space of the fade coefficient vector is a subspace of one of these 14 singular fade subspaces, i.e., 

the fade coefficient vector should belong to one of these 14 vector subspaces. For example, [2 1 + j]^ 

belongs to the vector subspace ( ) and is a deep fade matrix. 

\ 0.5 + 0.5J / 

Note that the singular fade subspaces depend only on the number of transmit antennas nt and the 
signal set S. They are independent of the number of receive antennas Ur, as illustrated in the following 
example. 

Example 2: Consider the 2 x 2 MIMO system with 4-PSK signal set S = {±l,±j}. The set of 14 
singular fade subspaces for this case is the same as that of 2 x 1 MISO system given in ([3]). For a 
fade coefficient matrix to be a deep fade matrix, both the rows should belong to one of these 14 vector 

1+i 1 r iT 



subspaces. For example. 



1 0.5 + 0.5j 
1 



is a deep fade matrix since 



2 1+i 



and 



1 0.5 + 0.5i 



belong to the vector subspace 

[0.5 + 0.5i 

The dimension of the singular fade subspace (Ax)-"-, and the transmit diversity order of the pair- 
wise error event (x — ;■ x'), are inherently connected, where Ax = x — x' and x, x' G 5"'. With spatial 
multiplexing, the transmit diversity order of the pair-wise error event (x — )• x') is 1 while di7Ti((Ax)-'-) = 
rit — l. It is the presence of these nj — 1 dimensional singular fade subspaces that results in a transmit 
diversity order of 1. 

The receive diversity order n,. comes due to the fact that for a fade coefficient matrix to be a deep fade 
matrix, all the rir rows of the fade coefficient matrix should belong to the same singular fade subspace. 

The use of full diversity space times space time codes results in the maximum transmit diversity order 
nt- In the next subsection, the connection between the singular fade subspaces of space time codes and 
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transmit diversity order will be established. 

B. Singular Fade Subspaces for the Collocated MIMO system with Space Time Coding 

Consider the case when Tx uses a space time code C of size nt x T, where T > iit. Let C(x) denote 
a codeword matrix of the space time code, where x € 5^, where K denotes the number of independent 
complex symbols transmitted. Similar to the spatial multiplexing case, the effective constellation at Rx 
which is a subset of C"'^^ can be defined. It is easy to verify that the minimum distance of the effective 
constellation at Rx becomes zero when Row{Yi) is a subspace of the vector subspace Col^ (C (Ax)) , 
for some Ax G A5^. Note that Co/-*- (C (Ax)) denotes the vector subspace {u : u'^v = 0,Vv S 
Col (C (Ax))}. The vector subspaces Col-^ (C (Ax)) are the singular fade subspaces for the nt transmit 
antenna system with the space time code C. 

Note 1: Even though the probability that RowCH) is a subspace of one of the singular fade subspaces is 
zero, with a non-zero probability i?ou;(H) falls in the neighbourhood of a subspace of one of the singular 
fade subspaces, which results in low values of the minimum distance of the effective constellation. 

The dimension of the singular fade subspace Col^ (C (Ax)) is equal to nt — rank(C (Ax)), while 
the transmit diversity order for the pair-wise error event (x — )• x'),x, x' G S^ equals rank(C(Ax)) 
lITSl . where Ax = x — x'. With every pair-wise error event (x — ^ x'), we can associate a singular 
fade subspace Col^ (C (Ax)) . Among all the pair-wise error events, those error events for which the 
codeword difference matrix has the least rank determine the overall system transmit diversity order. 
Equivalently, among all the pair-wise error events, those error events for which the associated singular 
fade subspace has the largest dimension will dominate the overall error probability. This is expected 
since among all the singular fade subspaces, the probability that Row(H.) falls in the neighbourhood of a 
subspace of the singular fade subspace, will be the largest for those singular fade subspaces which have 
the largest dimension. 

If the space time code is such that C (Ax) is full rank for all Ax 7^ Ok, all the singular fade 





subspaces Col (C (Ax)) collapse to be the zero-dimensional trivial singular fade subspace 
thereby ensuring a transmit diversity order of nt for all the pair-wise error events. 

Example 3: Consider the 2x1 MISO system with Alamouti space time code whose design matrix 
. Since the design matrix is full rank for all choices of xi and X2, the column 
difference matrix is always C^ and hence all the singular fade subspaces collapse 



is given by 

-X2 Xi 

space of the codeworc 



Xl X2 



to be the zero dimensional trivial singular fade subspace 



Equivalently, all the pair-wise error 
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events (xi,X2) —5- (x'-l,X2) have a transmit diversity order 2. The full rank Alamouti space-time code 
removed the effect of all the vector subspaces which were non-trivial singular fade subspaces for the 
spatial multiplexing system, thereby increasing the diversity order of all the pair-wise error events from 
1 to 2. 

Example 4: For a 2" x 2" Generalized Linear Complex Orthogonal Design (GCOD) |[T6l . the design 
matrix G2">(xi, X2, . . . Xq+i) constructed iteratively is given by, 

G2— i(xi,X2,. . .Xa) Xa+ll2— 1 

-x:+ll2a-i GH..(xi,X2,...xJ 

The codeword difference matrix for the GCOD is full rank for any signal set. Hence, irrespective of the 
signal set, the trivial singular fade subspace (02*) is the only singular fade subspace for the GCOD. 



Example 5: Consider the 4x4 Quasi-Orthogonal Design (QOD) lITTl . whose codeword matrix is given 



by 



xi 

X2 
X3 

X4 



-xz 



Xa 



Xi 



XT 
9 



X4 
-X3 

-X2 
Xl 



Let Ax, 



x[. Irrespective of the signal set used, the minimum rank 



of the codeword difference matrix for the 4x4 QOD is 2. For example, when Axi = Ax4 = Asi 
and Ax2 = — AX3 = As2, the rank of the codeword difference matrix is 2. Equivalently, there exists a 
non-trivial singular fade subspace, / Usi As2 -As2 Asi] , [-Asa Asl -Asl -Asa] ) • 

Note that the 2x2 Alamouti code removes the effect of the harmful non-trivial singular fade subspaces 
for any signal set. On the other hand, for the 4x4 QOD there exists non-trivial singular fade subspaces 
for any signal set. 

In general, a space time code can offer full transmit diversity for some but not all signal sets. In other 
words, for some signal set, a space time code might have only the trivial singular fade subspace, while 
for some other signal set, the same space time code might have non-trivial singular fade subspaces. For a 
space time code which does not offer full transmit diversity for a signal set, there would exist non-trivial 
singular fade subspaces. These are illustrated in the following example. 

Example 6: Consider the 2x2 Co-ordinate Interleaved Orthogonal Design (CIOD) |[T8l whose code- 
where xi,X2 G {±l,ibj}. Let Axj 



word matrices are of the form 



x^ + JX2 



xf + jx{ 



Xi 



The codeword difference matrix is not full rank in the following two cases: 

Case 1: Axf = Ax2 = and at least one out of Ax{ and Ax^^ is non-zero. For this case, the singular 
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fade subspace is given by 



Case 2: Ax{ = Ax^^ = and at least one out of Axf and Ax2 is non-zero. For this case, the singular 
fade subspace is given by 



and 



However, 



Hence, there exists the following two non-trivial singular fade subspaces: 

when the signal set is e-'^jitl, ±j}, where 9 is not a multiple of ^, the 2x2 CIOD offers full transmit 
diversity. Equivalently, there are no singular fade subspaces other than the trivial singular fade subspace 
for the 2 X 2 CIOD with the signal set e-'^{±l, itj}, when 6 is not a multiple of ^. 
Example 7: Consider the 4x4 CIOD ifTSl whose codeword matrices are of the form 



xf + jxl 


X2 -\- jx^ 








~^2 ~I"J^4 


xf - jxl 














xf + jx( 


xf + 








— X4 + JX2 


^3 ~ 



where xi,X2,X3,X4 G {ibl,ibj}. For the 4-PSK signal set considered, this STC does not offer full 
transmit diversity and there are pair- wise error events which have a transmit diversity order less than 2. 
The determinant of the codeword difference matrix for this S TC is given by. 



(|Aa;f I' + \Axif + \Ax^\^ + |A4|') (|Axf |' + |Aa:(|' + \Axf\^ + \Axi 

Hence the code-word difference matrix is not full rank in the following two cases: 
Case 1: Axf = Axg = Axf = AX4 = and at least one out of Ax(, Axf , Axf Axf is non-zero. 
For this case, the first two columns of the codeword difference matrices are zeros. The column span of 
the codeword difference matrix is ( [o 1 o] , [o 1 1 ) and hence the corresponding singular 
fade subspace is given by / h o] , [o 1 o] V 

Case 2: Ax{ = Axf = Ax2 = Axf = and at least one out of Ax{, AX3 , Axf Ax^ is non- 
zero. Similar to Case 1, it can be shown that the singular fade subspace for this case is given by 

( [0 1 0] , [0 1] y 

Hence, for the 4x4 CIOD, with 4-PSK signal set xi,X2,X3,X4 G {±l,ibj}, there exists two non- 
trivial singular singular fade subspaces. Similar to the 2x2 CIOD, when the signal set is a rotated 4-PSK 
signal set, €•'^{±1, itj}, where 6 is not a multiple of |, the 4x4 CIOD offers full transmit diversity 
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and there are no non-trivial singular fade subspaces. 

In general, the 2x2 CIOD given in Example [6] and the 4x4 CIOD given in Example |7J offer full 
diversity for those signal sets for which the Co-ordinate Product Distance (CPD) \j is non-zero lITSl . 
Equivalently, there are non-trivial singular fade subspaces for the 2x2 and 4x4 CIOD, for signal sets 
whose CPD is non-zero. In fact, this is true for any Generalized Co-ordinate Interleaved Orthogonal 
Design (GCIOD), as illustrated in the next example. 



|[T8l whose codeword design matrix is given by. 



The complex 



Example 8: Consider the 2" x 2" Generalized Co-ordinate Interleaved Orthogonal Design (GCIOD) 

G2a-l{xi,. .. ,Xa) 

6*20-1 (ia + l, . . . ,X2a) 

number Xj = x^ + jxf ^ , where (r)^ denotes r modulo s and G20-1 (xi, . . . , Xa) is the codeword 
matrix of the GCOD lfT6l of size 2""^. The determinant of the codeword difference matrix is given by 
(Er=i(l^^?l' + \^^la+^)2a\'') (E"=i(|Aa:fr + \^^fa+^)2S) ■ Thc detcmunant is non-zero for those signal sets 
for which the CPD is non-zero and there are no non-trivial singular fade subspaces. For those signal sets 
for which the CPD is zero, the determinant becomes zero under the following two cases: 

Case 1: Ax^ = Ax^ .^^ =0,Vi G {!,••• , a} and at least one of the elements of the set {Aa;f,Aa;[^+j)2^, 
1 < i < a} is non-zero. It can be verified that the singular fade subspace for this case is given by 
(ei,e2,... ,620-1) , where Cj denotes the 2" length vector whose i^^ component is one and all other 
components are zeros. 

Case2: Axf = Axf , j-, = 0, Vi G {1, . . . , a} and at least one of the elements of the set {Aa;!, Aa;^^^,)^^, 
1 < i < a} is non-zero. For this case, the singular fade subspace is given by (620-14-1, • • • , 62") . 

III. Singular fade subspaces for the two-way relaying scenario 

In the previous subsection, the notion of singular fade subspaces was introduced and its connection 
to the transmit diversity order of the MIMO system with collocated antennas was established. Since the 
MA phase of the two-way relaying scenario can be viewed as a virtual 2x1 MISO system, there exists 
singular fade subspaces for this case as well. 

In Subsection IIII-A[ the singular fade subspaces for the two-way relaying scenario are identified. The 
reason why the adaptive network coding schemes based on the DNF protocol proposed in |J6] and lITTI - 
|[T2l mitigate the effect of these harmful singular fade subspaces is discussed. In Subsection IIII-B[ it is 

'The CPD between two complex numbers x and y is defined to be \x^ — y^\\x' ^ y'\- The CPD of a signal set is defined 
to minimum among all CPDs between pairs of points in the signal set 1181 . 
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shown that minimizing the harmful effect of these singular fade subspaces can also be achieved by a 
proper choice of the DSTC, without any need to adaptively change the network code at R according to 
channel conditions. 



A. Singular Fade Subspaces for the DNF Protocol 

Let AxA = XA — x'j^ and Aj;^ = xb — x'g G A5. From the discussion in Section JIl it follows that 

The ratio 



Asi 



1 

-Ax A 



the singular fade subspaces for the DNF protocol are of the form 

~^^^ determines all the singular fade subspaces for the DNF protocol. In lITTI - lITll . the ratio ^^^ 
was called the singular fade state. 



As mentioned earlier in Section |II1 dim 



1 

-Ax/ 



and the diversity order for the pair-wise error 
event that a pair {xa.xb) is wrongly decoded at R as (x^,x'^) (denoted as {xa^xb) -^ {x'^.x'^)) 

are inherently connected. The diversity order for the error event {xa,xb) 

1 



{x'j<^,x'q) is equal to 

/ / 1 \ \ 
rank([Axyi Ax^]) = 1 while dim 

\ \ [ AXB J ' / 

Let Sji{hA,hB) = {Haxa + It-bxb '■ xa,xb G S} denote the effective constellation at R. Let 
d-minihA, ^b) denote the minimum distance of Sji{hA, Hb)- When [hA /i_b]^ falls in one of the singular 
fade subspaces, dmin{hA, h-B) becomes zero. Even though the probability that the vector [/ia/ib]"^ belongs 
to a singular fade subspace is zero, dmin{hA, Iib) is greatly reduced when [Iia Iib]^ falls close to a 
singular fade subspace, a phenomenon referred as distance shortening. For /S.xa t^ and /S.xb t^ 0, 
the CNC algorithm |[6l avoids the distance shortening occurring in the neighbourhood of a singular fade 



subspace 



, by ensuring that ^^'^'^^{xa,xb) = fJ,^-'^'^"{x'j^,x'^), i.e., R does not distinguish 
the pairs {xa,xb) and {x'y^,x'^), which are said to be clustered together In fact, for every realization of 
[liAhB] (not necessarily in the neighbourhood of singular fade subspaces), the CNC algorithm chooses 
the network coding map which results in the best distance profile at R by appropriate clustering of the 
signal points. The scheme proposed in lITTI - |[T2l avoids distance shortening in the neighbourhood of 
singular fade subspaces by proper choice of clustering for only the singular fade subspaces and not for 
every realization of the channel fade coefficients. 



Consider the two singular fade subspaces: 



AxB 



and 



AxA 





The 



distance shortening which occurs in the neighbourhood of these singular fade subspaces is unavoidable, 
since the pairs {xa,xb) and {xa,x'^) (and also the pairs {xa,xb) and {x'y^,XB)) which result in these 
singular fade subspaces cannot be clustered together without violating the exclusive law. Such singular 
fade subspaces are referred as the non-removable singular fade subspaces. The dimension of these singular 
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fade subspaces is one or equivalently, the error events {xa,xb) — )• {xa,x'q) and {xa,xb) -^ {x'^,xb) 
always result in diversity order one. The singular fade subspaces other than the non-removable singular 



fade subspaces are referred as the removable singular fade subspaces. 

1 



The removable singular fade subspaces are of the form 



-Aja 



,AxA / 0,Axs / 0, which are 



and 
e singular fade 



dependent on the signal set S used. The non-removable singular fade subspaces are 
which are independent of the signal set used. Owing to the presence of non-removab 
subspaces, the overall diversity order of the DNF protocol cannot exceed one. 

From the discussion above, it is clear that there are two classes of singular fade subspaces: removable 
and non-removable. The non-removable singular fade spaces are created by the channel and is independent 
of the signal set used. Whatever may be the choice of the network code, the harmful effects of these 
non-removable singular fade subspaces cannot be mitigated. The harmful effect of the removable singular 
fade subspaces, which are created by the signal set, can be removed by a proper choice of the adaptive 
network coding map at R, as in 161 and ifTTI - |[T2l . 

To sum up, in the DNF protocol, the transmissions from the nodes A and B are allowed to interfere at 
R and the effect of MAI is effectively mitigated by adaptively changing the network coding map, thereby 
removing the harmful effect of all the removable singular fade subspaces. 



B. Singular Fade Subspaces for the DSTC Scheme 

Let Axa = xa — x'^ and Axb = xb — Xg G A5^. Then C(Axa, Axb) = C(xa,xb) — C(x'^,Xg) 
denotes a codeword difference matrix of the DSTC, where C(xa,xb) is the codeword matrix of the 
DSTC defined in ([T]). From the discussion in Section JIJ it follows that the singular fade spaces for the 
proposed DSTC scheme are of the form Col^ (C (Axa, Axb)) • 

Consider the singular fade subspaces of the form Co/-*- (C (O2, Axb)) and Co/-*- (C (Axa, O2)) , 
where Axa, Axb 7^ O2. The first row of the matrix C(02, Axb) has both the entries to be zero. 



Hence, Co/ (C (02,Axb)) 



and the singular fade subspace Col-^ (C (O2, Axb)) 



By a similar reasoning, Col-^ (C (Axa, O2)) = ( 
If the DSTC codeword matrices are such that rank 



(C(AxA, Axb)) = 2, VAxa / O2 and Axb / O2, 



all the singular fade subspaces Col^ (C (Axa, Axb)) collapse to be the trivial singular fade subspace 
(O2) . Equivalently, all the pair- wise error events C(xa,xb) — )■ C(x'^,Xg),XA / x'^,xb 7^ Xg, have 
diversity order 2. Hence, for a properly chosen DSTC, other than the trivial singular fade subspace, the 
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singular fade subspaces are only the two non-removable singular fade subspaces, while for the DNF 
protocol, in addition, we had the removable singular fade subspaces. In this way, by a proper choice 
of DSTC, the occurrence of the removable singular fade subspaces is avoided at the transmitting nodes 
itself, without any CSIT. 

Hence, we have the following design criterion referred as the singularity minimization criterion for 
DSTCs for two-way relaying: The DSTC codeword difference matrices C(Axa, Axb) need to be full 
rank for all Axa 7^ O2 and Axb 7^ O2, to minimize the number of singular fade subspaces. DSTCs 
satisfying the above criterion are referred as the singularity minimal DSTCs. 

Hence for a DSTC which is singularity minimal, the only error events which result in diversity order 1 
are of the form C(xa,xb) — ?> C(xA,Xg),Xg / xb and C(xa,xb) — )■ C(x'^,xb),x'^ / xa- Hence, 
the overall coding gain is equal to minimum among all the non-zero singular values of the codeword 
difference matrices which are of the form C(02, Axb) and C(Axa,02) llT5l . Note that the matrices 
C(02, Axb) and C(Axa,02) are of rank 1 and have only one non-zero singular value. We have 
the following coding gain criteria for singularity minimal DSTCs: the minimum among all the non-zero 
singular values of the codeword difference matrices which are of the form C(02, Axb) and C(Axa, O2) 
needs to be maximized. 

Example 9: Consider the DSTC . This DSTC is nothing but the scheme where A and B 

XB^ 

transmit in separate time slots, making sure that their transmissions do not interfere at the relay. Even 
though this DSTC avoids all the removable singular fade subspaces, the end-to-end rate in complex 
symbols per channel use is less than that of the DNF protocol. 

C. A Construction of Singularity Minimal DSTCs for Algebraic Signal Sets 

A signal set is said to be algebraic if all the signal points of the signal set are algebraic numbers over 



i 



All the commonly used signal sets like QAM and PSK are algebraic signal sets. In this subsection, 

a b 



be 

c d 

the form 



a class of DSTCs which are singularity minimal for algebraic signal sets is provided. Let 

a full rank complex matrix. Consider the class of DSTCs whose codeword matrices are of 

_,. , a{xAi+e^XA2) b{xAi + e^XA2) 

C(xa,xb) = 

c{xbi +e^XB2) d{xBi +e^XB2) 
Proposition I: The class of DSTCs whose codeword design matrices are of the form given above are 

singularity minimal for all algebraic signal sets. 

^A number is said to be algebraic over Q if there exists a polynomial with coefficients from Q of which the number is a 
root. If there does not exist a polynomial with coefficients from Q of which the number is a root, the number is said to be 
transcendental 1 191 . 
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Proof: The proof is as follows: For Axa / O2 and Axb / O2, at least one of the two components of 
Axa as well as Axr should be non-zero. Hence, (Ax^^+e-^'Ax^J / and (Axb^+c^Axb^) 7^ 0, since 
e^ is transcendental G whereas Ax^^, Ax^^, Ax^j and Ax^a are algebraic over Q. The codeword differ- 
ence matrix C(Axa, Axb) is full rank for all Axa / and Axb 7^ 0, since det(C(AxA, Axb)) = 



{ad - 6c)(AxAi + eJAxAj(A2;B^ + e^AxsJ / 0. 



Example 10: Consider the case when 



a b 
c d 



1 
1 



Let 4-PSK be the signal set used at A and 



B. The DSTC codeword matrix for this case is given by, C(xa,xb) 



{xAi + e^XA^) 

{xBi + e^XB2) 

and B are made to transmit in two different time slots which results iii low decoding complexity at R, 



since A's and B's transmissions can be decoded independently. It can be verified that the coding gain for 
this DSTC is approximately 0.6877. 



Example 11: Consider the case when 



a h 
c d 



1 1 
-1 1 



Let 4-PSK be the signal set used at A and B. 



-(a;si + e^XB-i) {xbi + e^ssj) 



The DSTC codeword matrix for this case is given by, C(xa, xb) = 4= 

The scaling factor of 4= is to ensure unit average energy per symbol per time slot. It can be verified that 
the coding gain for this DSTC is approximately 0.6877, same as that of the DSTC given in Example [TOl 
The coding gain of the DSTCs given in Examples [TOl and [TT] is approximately 0.6877, which is less 
than the minimum distance of the unit energy 4-PSK signal set, which is \/2. In the next section, it is 
shown that for DSTCs over square QAM and 2'*'-PSK signal sets, the coding gain is upper bounded by 
the minimum distance of the signal set and explicit DSTC constructions which achieve this bound with 
equality are provided. 



IV. Singularity Minimal, Coding Gain Maximal DSTCs over QAM and PSK signal sets 

In this section, it is shown that the coding gain of the DSTCs over square QAM and 2^ -PSK signal 
sets are upper-bounded by the minimum distance of the signal set. In Subsection IIV-A[ a condition under 
which a singularity minimal DSTC over square QAM and 2'^ -PSK signal set meets the upper bound with 
equality is obtained and explicit constructions of DSTCs are provided. In Subsection |IV-B[ the constructed 
DSTCs are shown to be fast ML decodable, i.e., the ML decoding complexity of the constructed DSTCs 
is shown to be less than the brute-force decoding complexity which is 0{M^). 

Note that the generator matrices Ma and Mb at A and B should be such that the average energy per 
time slot is unity, i.e., E(|| xaMa f ) < 2 and E(|| xbMb P) < 2. 

'By Lindemann-Weierstrass theorem 1191 , e-" is transcendental for all g G Q. 
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Lemma 1: For singularity minimal DSTC over S, where 5 is a square QAM or 2^-PSK signal set, 
the coding gain is upper bounded by the minimum distance of the signal set S. 

Proof: See Appendix A. ■ 

In the following subsection, the condition under which the upper-bound given in the previous lemma is 
satisfied with equality is identified and explicit construction of DSTCs are provided. 

A. Constructions of Singularity Minimal, Coding Gain Maximal DSTCs over QAM and PSK signal sets 

The following proposition states that for DSTCs over S, choosing Ma and Mb to be unitary matrices 
ensures that the upper-bound on the coding gain is satisfied with equality, for QAM and PSK signal sets. 

Proposition 2: For singularity minimal DSTCs over square QAM or 2'*'-PSK signal sets, the coding 
gain is maximized when the generator matrices Ma and Mb at A and B are unitary matrices. 

Proof: When Ma and Mb are unitary matrices, || AxaMa || = || Axa || and also || AxbMb || = 
II Axb II • Hence, min || AxaMa ||= min l^^;^, I = dmmiS) and similarly 

AxAeA52, AxAiSA^.AxAi^O 

AXA7^02 

min II AxbMb II = dmin(S), where dminiS) denotes the minimum distance of S. 

AxBGA5^ 
AxBT^Oa 

The coding gain of the DSTC is the minimum among all the non-zero singular values of the codeword 
difference matrices which are of the form C(02, Axb) and C(Axa, O2), i.e., the coding gain is equal 

to min < min || AxaMa ||, min || AxbMb || > , which is equal to dmin{S)- ■ 

AxaSAS^, AxoeAS^, 

l^ Axa 7^02 Axb 7^02 ) 

In the following examples, constructions of singularity minimal DSTCs whose generator matrices are 



Construction 1: Consider the DSTC over S for which Ma = A= and Mb = A= 





- 


- 


1 


a. 


a 


^5 


a(j> 


a(j> 



unitary are provided. 

ja a 
jctcj) caf) 

where = ^~^2 ' ^ ~ ^~2 ' ^ — ^ ~^ J ~ 3't> and a = 1 + j — j0. The DSTC codeword matrix is of 

the form C(xa,xb) = . The codeword difference matrix C( Axa, Axb) is full rank for all 

xbMb 

Axa 7^ and Axb 7^ 0, when the signal points belong to Z[j] |[20l . Hence the DSTC is singularity 
minimal for all signal sets whose signal points belong to Z[j]. Also, since Ma and Mb are unitary, for 
square QAM signal set, the DSTC maximizes the coding gain. 

Note 2: The DSTC given in Construction [T] was constructed in EOl towards satisfying the design 
criterion formulated in ^T\\ for the two-user non-cooperative Multiple Access Channel (MAC). In f20l . 
the DSTC given in the above example was shown to be DMT optimal for two-user MAC. 
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where 



Construction 2: Consider the DSTC for which Ma = I2 and Mb = 

(pg = tan~^ 1/5. The DSTC codeword matrix C(xa,xb) is given by, 

XAi XA2 

xbi cos <j>g + XB2 sin 09 e-* (— x_Bi sin cj>g + XB2 cos ( 

For a complex number a, let Q(a) denote the smallest field containing Q and a. It is shown in Lemma 
|2] below that choosing ^ = f ensures singularity minimality for signal sets (for example QAM) whose 
signal points belong to Q(j) and choosing = -^ ensures singularity minimality for signal sets (for 
example 2^-PSK) whose signal points belong to Q{e-'^). Also, since Ma and Mb are unitary, this 
DSTC maximizes the coding gain, for square QAM and 2^-PSK signal sets. The advantage of this 
construction over Construction [T] is that encoding at node A is simple, since it does not involve any 
linear combination of xai and XA2 ■ 

Lemma 2: For the DSTC given in construction ID choosing 9 = j ensures singularity minimality for 
signal sets whose points belong to Q(j) and choosing 9 = ^ ensures singularity minimaUty for signal 
sets whose signal points belong to <^{e'^). 

Proof: The proof is given for the case when the signal points belong to Q(j). The proof for the 
case when the signal points belong to ^{e^^) is exactly similar and is omitted. 

Let /S.XA, = XAi ~ ^'Ai ^^^ ^^B, = ^B, — x'^., where XAi,x'^,,XB,,x'^^ G 5 C Q(j). and i G {1,2}. 
To prove singularity minimality, it needs to be shown that when at least one out of Axai ^^id Ax^^ 
(Axb, and Ax^,) is non-zero, the codeword difference matrix is full rank. The ratios . "°^ and — . "°^ 
belong to Q(j) while tan^g = \/5 does not belong to Q(j). Hence, Axb^ cos (pg + Axb^ sin (pg / 
and —AxBiSm(pg + Axb2 cos (pg 7^ 0. Since sm(f)g = -^ and cos(t)g = -y=, Axa2{Axbi cos (pg + 
AxB2sin(l)g) and AxAi{—AxBiSin(pg + AxB2Cos(l)g) belong to Q(j, a/S, \/6), where Q(j, \/5, -v/6) 
denotes the smallest filed containing Q, j, \/5 and ^/6. The determinant of the codeword difference 
matrix is given by, 

AxAie^~{-AxBi sm(f)g + Ax^^ cos(pg) - Axa2{Axbi cos(pg + Axb^ sin^g). 

The determinant is non-zero since the ratio ^^2^ xg^cosipg xg^sm^gj belongs to Q(i, 1/5, -\/6), while 
e-' 4 does not belong to Q(j, -\/5, a/6). ■ 
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B. Decoding Complexity of Singularity Minimal, Maximal Coding Gain DSTCs over S 

After the two MA phases, R jointly decodes for the two message vectors xa and xb of A and B 
respectively. In general, the complexity of this joint ML decoding at R is 0{M'^), where M is the 
cardinality of the signal set S. The choice of the generator matrices Ma and Mb being unitary not 
only maximizes the coding gain for QAM and PSK signal sets, but also results in a reduced decoding 
complexity at R. 

The following proposition states that when conditional ML decoding ll22l . ll23l is employed, the 
decoding complexity of the DSTCs constructed in the previous section for which the generator matrices 
Ma and Mb are unitary is 0{M^) for any arbitrary signal set and is 0{]VP) for square QAM signal 
set. Note that the brute force decoding complexity is 0{M^). 

Proposition 3: When the generator matrices of the singularity minimal DSTC over S are unitary, the 
decoding complexity using conditional ML decoding is 0{M^) when the signal set S is arbitrary and is 
0{M^) when the signal set S is square QAM. 

Proof: See Appendix B. ■ 

Compared with the DNF protocol, the decoding complexity is more for singularity minimal coding 
gain maximal DSTCs over S. For the DNF protocol, the decoding complexity is 0{M'^) for non square 
QAM signal sets while it is 0{M) for square QAM signal setQ As indicated by the simulation results in 
the next section, the proposed DSTC offers slightly better performance than the adaptive network coding 
scheme and eliminates the need for adaptive switching of network coding maps at R. But this comes at 
the cost of increased decoding complexity at R. 

V. Simulation Results 

All the simulation results presented are for the case when the end nodes use 4-PSK signal set. By 'DSTC 
1' and 'DSTC 2' we refer to the DSTCs given in Construction [T] and Construction |2] respectively. As a 
reference scheme, we consider the scheme in which XOR network code is used irrespective of channel 
conditions and no DSTC is employed, which is referred as 'XOR NAV code'. Assuming unit noise 
variances at all the nodes, the average energies of the transmissions at the nodes, which are assumed to 
be equal, is defined to be the Signal to Noise Ratio (SNR). The proposed DSTC scheme is also compared 
with the adaptive network coding schemes proposed in ||6] and lITTI - |[T2l . Since for 4-PSK signal set, the 
adaptive network coding scheme based on the Nearest Neighbour Clustering (NNC) algorithm proposed 

■*For the DNF protocol, with QAM signal set, conditioning on xa, xb can be decoded with constant decoding complexity 
by rounding off to the nearest integer, which results in an overall decoding complexity of 0{M). 
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in ||6l and the scheme based on Latin Squares proposed in lITTI - |[T2l turn out to be the same, without 
distinguishing them we refer to both as 'adaptive NAV code'. Fig. |2] shows the SNR vs BER performance 
for different schemes for the case when all the fading coefficients are i.i.d. and Rayleigh distributed. In 
Fig. [3] and Fig. |4] similar plots are shown for a Rician fading scenario with Rician factors |j of dB and 
5 dB respectively. From Fig. |2]|4l it can be seen that the diversity order is one for all the schemes. Also, 
it can be seen that at high SNR, both 'DSTC 1' as well as 'DSTC 2' offer nearly the same performance 
and they perform better than the 'XOR NAV code' as well as the 'adaptive N/W code'. For a Rayleigh 
fading scenario, at high SNR, the DSTCs offer a gain of 2 dB over 'XOR NAV code' while the 'adaptive 
NAV code' offers a gain of about 0.5 dB over the 'XOR NAV code'. For a Rician factor of dB, at high 
SNR, the DSTCs offer a gain of 2 dB over 'XOR N/W code' while the 'adaptive N/W code' offers a gain 
of about 1.2 dB over the 'XOR N/W code'. For a Rician factor of 5 dB, at high SNR, the DSTCs offer 
a gain of 5.5 dB over 'XOR N/W code' while the 'adaptive N/W code' offers a gain of about 4 dB over 
the 'XOR N/W code'. The reason why the DSTC based scheme performs better than the adaptive N/W 
coding scheme is as follows: during the BC phase always a 4 point signal set is used for the DSTC based 
scheme, while depending on channel conditions 4 point or 5 point signal set is used for the adaptive 
network coding scheme ||6], lITTI . 

VI. Discussion 

A DSTC scheme was proposed for the two-way relaying scenario. It was shown that deep channel 
fades occur when the channel fade coefficient vector falls in a finite number of vector subspaces called the 
singular fade subspaces. The connection between the dimension of these vector subspaces and the transmit 
diversity order was established. Design criterion to minimize the number of singular fade subspaces for 
the DSTC scheme and maximize the coding gain were obtained. Explicit low decoding complexity 
constructions of DSTCs were provided. The problem of constructing singularity minimal DSTCs with 
decoding complexity same as that of the DNF protocol, without sacrificing the coding gain, remains 
open. Extending the DSTC scheme for two-way relaying with multiple antennas and multi-way relaying 
are possible directions for future work. 
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Appendix A 
Proof of Lemma 1 

Since IMaM^ is Hermitian, it is unitaiily diagonalizable, i.e., MaM^ = UaAaU^, where Ua is a 



unitary matrix and Aa is a diagonal matrix with diagonal entries denoted as A^^ and Xa^ ■ ^ 

an ai2 



and A^, are non-negative since IvIaM^ is positive semi-definite. Let IMa 



021 022 



ote that Aai 
We have, 



2,1 |2 

a2i\ + a22 , 



E(|| xaMa f) = \aii\'E{\xAA) + \ai2rH\^A2n + \a2i\M\^AA) + \a22m\xA,n = la 

since E{xa^x\J = E{xa2x\J = for square QAM and 2^-PSK signal sets. Since E(|| xaMa |P) < 2, 

we have janp + |ai2p + |a2ip + |a22p = Trace(MAMA^) = XA^ + Aa, < 2. 

The coding gain of the DSTC is the minimum among all the non-zero singular values of the codeword 
difference matrices which are of the form C(02, Axb) and C(Axa,02), i.e., the coding gain is equal 

to min < min || AxaMa ||, niin || AxbMb 

I AxaSAS^, AxbGAS^. 

l_ AxAT^Oa Axb 7^02 

Let dmin{S) denote the minimum distance of the signal set S. 

Consider || AxaMa P= AxaMaMa^Axa^^ = AxaAaAxa^ = \a^\^xa,? + A^JAxaJ^, 
where Axa = AxaUa — [Ax^i Ax^J- 

Let UAi = [uAii UA12] and UA2 = [UA21 ua-^-i] denote the rows of Ua- For Axa = [^xa^ 0], 
II AxaMa f= \/^xa,?{\ua,^?\a^ + I^AiJ^AaJ. 

Hence, we have, min || AxaMa ||^< (frnin{S){\uAii\^ \ai + \uai2?\a2)- Similarly, we have, 

AxaSAS^, 
Axa/02 

min II AxaMa ||^< d^i„(5)(|uA,i I^Aai + |ua„ I^Aaj). Since Ua is unitary 1^^^ P = I^a,, P and |uai2 P 

AxaSAS^, 
AXA#02 

|t(A2iP- Therefore, we have, 

min II AxaMa 11^ < rfLn(5)min{(|ttAii|^AAi + |«Ai2l^AA2),(hAii|^AA2 + ^Aial^AAi)}. (4) 

AxaSAS^, 
Axa ^02 

Since Ua is unitary, |mAiiP = 1 — I^Ai2p- For ^ given Aai and Aa2, the upper-bound in (|4]) is 
maximized over all IuauP when the two terms inside min are equal, i.e., IuaiJ^Aai + |^Ai2pAA2 = 
KaiiPAa2 + KAi2pAAi, for which Iua^P = \ and this maximum value is equal to d^^^{S) '^^ ^ "^^ . 
Since, Aai + Aa, < 2, the maximum value of the upper-bound in dUl is less than or equal to (i,^^^^„(5). 
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Hence, min AxaMa < dminiS). Similarly, it can be shown that min AxbMb is 

AxBeA5^ AxAeA5^ 

Axb^02 Axb^02 

also upper-bounded by dmin{S). Hence, the coding gain of the DSTC over square QAM or 2'^-PSK 
signal set is upper-bounded by dmin{S)- This completes the proof. 



Appendix B 
Proof of Proposition 3 

To prove the proposition, we adopt a procedure similar to the one used in |[23l . 

Let YR = [y^^ y^^^ y^^ yj^f, St = [x^; x^^^ x^^ x;^^ xg^ x^^ xf ^ x;^J^ and ia = [z^^ 4^ z^^_ zj^f. 
The vector yR can be written as y^ = HgqX + zr, where Hgq is a 4 x 8 real matrix whose entries are 
functions of Ha and hs, determined by the DSTC. Using QR decomposition, the matrix Heq can be 
decomposed as Hgq = QR, where Q G M^^"^ is a orthogonal matrix and R € M"^^® can be written as 
[Ri R2], with Ri, R2 € M^^^, Ri being an upper-triangular matrix. The joint ML decoding metric at 
R is given by || yR - HgqX || = || Q'^yR - Rx || = || yR - Rx ||, where y'^ = Q'^yR. 

For a singularity minimal DSTC over 5, let the generator matrices be Ma = Ua and Mr = Ur, 
where Ua and Ur are unitary matrices. Let uai and ur; denote the i^^ rows of Ua and Ur 
respectively. Then the weight matrices of the DSTC defined in Q are given by, W^ = jW^ = 



UAi 



and Wg. 



jWl 



ilarly, W^^W^^ 



R H 



Bi — 

-j 





0^ 



UB: 



We have, W^^W^^ 



I H 



UAi 




JuS, O2 


= 


J 




and sim- 



Tl H 



R H 



. Hence, W^^W^^ + W^^W^^ = O2, where O2 denotes the 



2x2 null matrix. Also, W^ W^ = O2, since uai and UAa are orthogonal vectors. Hence, 
W^ W^ + W^ W^ = O2. Similarly, using the fact that Ua and Ur are unitary matrices, 
it can be shown that the following pairs of matrices are also Hurwitz-Radon orthogonaO {W^ , W^ }, 
{wi^, w^j, {wi^, wij, {w^,, wij, {Wg^, W^J, {Wg^, Wgj, {W%^, W^J, {W^^, 
WgJ, {W^^,W^J, {W|^,W^J. 

Let Xj denote the i*^ component of the vector x. The i*'^ and j*'^ columns of Hgq are orthogonal and 
hence the (i, j)*'^ entry of R (i < j) is zero for all realizations of /i^ and /i^, if and only if the weight 
matrices of the DSTC corresponding to the symbols Xj and Xj are Hurwitz-Radon orthogonal (follows 
from Theorem 2, ll23l ]l Hence the matrix R is of the form given below. 



•^Two matrices Mi and M2 are said to be Hurwitz-Radon orthogonal if MiM" + M2Mf = 0. 

^Theorem 2 in 1231 proves only the 'if part. However, following an approach similar to the proof given in 1231 , it is easy to 
show that the weight matrices of the DSTC corresponding to the symbols Xi and Xj need to be Hurwitz-Radon orthogonal, for 
the {ijjY'' entry of R (i < j) to be zero for all realizations of Ha and hs, and hence the 'only if part also holds. 
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R 



* 
0*00 
0*0 
0* 



(5) 



w^w^^+wgwi^" 



UAi ug^ 



UBiU^^ 



Note that * denotes possible non-zero entries. The claim is that all the entries denoted by * are non- 
zeros. It is clear that all the diagonal entries are non-zeros. For the (1,5)*^ entry in ([5]) to be a zero, 

0, which implies that uai and ub^ are orthogonal 

vectors. Then the vector ubi should belong to the one-dimensional subspace which is orthogonal to uai • 

Since UAa also belongs to this one-dimensional subspace and both ub^ as well as UAa are of unit norm, 

ubi = e-'^UAa, for some angle 9. In that case, the DSTC codeword difference matrix is of the form 

^ ^ , which is not full rank when A.XA2, ^xbi 7^ 0, Ax^^^ = Ax^^ = and hence 

the singularity minimization criterion is violated. Hence, (1, 5)*'^ entry shown by * in ([5]) is non-zero. By 
a similar argument, it can be shown that the other non-diagonal entries denoted by * in ([5]) are non-zeros. 

From the matrix R given in (l5]l, it can be seen that conditioning on the variables xb^ and XB21 the 
symbols xai and XA2 can be decoded independently |[23|| . Since the total number of choices for xbi and 
XB2 is M^ and independently decoding xa^ and x^^ requires 2M computations, the decoding involves 
2M^ computations and hence the decoding complexity at R is 0{M^). 

For square QAM signal sets, the decoding complexity can be further reduced, since the real and 
imaginary parts independently take values. From (jSj, it can be seen that conditioning on xb^ and x^^ , 
the real and imaginary parts of xa^ as well as xa^, can be decoded independently. Since decoding the 
real and imaginary points of a signal point in QAM signal set is of constant complexity independent of 
M(decoding can be done by rounding off to the nearest integer ||23l ). the ML decoding complexity is 
0{M'^) for square QAM signal sets. This completes the proof. 



February 26, 2013 



DRAFT 



26 



h 



h 



B 
x_ 



h' 



h' 






(a) MA Phase 
Fig. 1. Wireless two-way relaying 
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Fig. 2. SNR vs BER for different schemes for 4-PSK signal set for a Rayleigh fading scenario. 
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Fig. 3. SNR vs BER for different scliemes for 4-PSK signal set for a Rician fading scenario witfi a Rician factor dB. 
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Fig. 4. SNR vs BER for different scliemes for 4-PSK signal set for a Rician fading scenario with a Rician factor 5 dB. 
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